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Ordered phases on curved substrates experience a complex interplay of ordering and intrinsic 
curvature, commonly producing frustration and singularities. This is an especially important issue in 
crystals as ever-smaller scale materials are grown on real surfaces; eventually, surface imperfections 
are on the same scale as the lattice constant. Here, we gain insights into this general problem 
by studying two-dimensional smectic order on substrates with highly localized intrinsic curvature, 
constructed from cones and their intersections with planes. In doing so we take advantage of fully 
tractable "paper and tape" constructions, allowing us to understand, in detail, the induced cusps 
and singularities. 

PACS numbers: 61.30.Jf, 61.30. Dk, ll.lO.Lm 



I. INTRODUCTION AND SUMMARY 

Dramatic progress has been made in algorithmic 
origami; it is now possible to design nearly arbitrary 
three-dimensional constructions out of unstretched [1 or 
nearly unstretched [2 plaquettes, isometric to pieces of 
the Euclidean two-plane. At the same time, there has 
been theoretical and technological interest in crystalline 
and liquid crystalline order on curved substrates [3]-[7]. 
The latter problem can be studied in reduced complex- 
ity by considering surfaces with vanishing Gaussian cur- 
vature except at isolated points and curves. Were we 
to consider only intrinsic interactions between the sub- 
strate and the ordered phase, it follows that the in-plane 
positional and orientational order would be completely 
determined by the folds and conical points of the sub- 
strate. The connection between smectic textures, geo- 
metric optics m [8], and shocks |9] on flat and curved 
surfaces makes studies of the smectic phase amenable to 
exact analytic study [101 HI] while capturing the salient 
features of broken translational and rotational invariance 
[12 . Combined with the simplified geometries we con- 
sider, we are led to highly tractable models of order on 
curved backgrounds. 

Here we study equally-spaced smectic textures on in- 
finite cones, cones interesecting with planes, and cones 
intersecting with cones, the latter two standing in for 
simple bumps on surfaces and saddle-like regions, respec- 
tively; see Fig. [l] Not only do we see the occurrence of 
focal lines and cusps in the ensuing smectic textures, but 
we also see violations of the rules that usually govern 
the schlieren textures of the sample. Recall that curved 
geometry interacts with topological defects [13^ to alter 
the conservation of topological charge in much the same 
way the background intrinsic curvature changes the rules 
for the sum of the interior angles of a triangle. In the 
case of schlieren textures in flat geometries, defects mark 
the confluence of an even number of dark brushes. How- 
ever, even this simple counting rule is violated on curved 
surfaces, as we will demonstrate. 

We commence with a cone C embedded in R^. C is 




FIG. 1: Substrates with conical defects. Notice that all of 
these surfaces are intrinsically flat almost everywhere - they 
fail to be flat at points and lines only. 



a singular surface which is flat everywhere except at its 
apex, where all the Gaussian curvature is concentrated. 
The geometry of C may be conveniently examined by 
cutting the cone along a radial line L and laying it flat 
on a plane (Fig.|2|. This way, C looks like a disk with 
a circular sector of angle 5 removed and with its two 
straight edges identified. The angle 5 is called the deficit 
angle. A direct application of the Gauss-Bonnet Theorem 
shows that 5 is also the total Gaussian curvature of any 
region of the cone which contains its apex. 

In order to establish notation, let X, Y, Z and x^y he 
Cartesian coordinates as in Fig. [2j A parametrization of 
C is given by: 

X{R,^) =i?cos$sinC, 
Y{R,^) =i?sin<I)sinC, 
Z{R,^) =-RcosC, 

where, as can be easily seen, the apex angle 2( is related 
to S by sinC = ^f^. In terms of x^y (see Fig. [2]) and 
their polar coordinates, r = y^x^^^y^ and tanc/) = y/x^ 
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FIG. 2: A cone is isometric to a cut planar disk with two 
radial lines (denoted by L above) identified. The point A is 
the apex of the cone. 



we have 



R = r, 
^ = 



(1) 



II. BUILDING THE LAYERS 

A necessary condition for the layers to be uniformly 
spaced is that their normal vector field points along 
geodesies of the surface |8l[14]. In the flattened model, 
these geodesies are just the straight lines of the plane. 
We know that defects (even in flat space) tend to con- 
centrate on lower dimensional sets in order to save energy 
so, in a 2-dimensional surface, this means that point de- 



fects are favored and this gives rise to layered structures 
in the form of wavefronts emanating from a point. Note 
that the case where this point is taken to infinity for- 
mally corresponds to a defect free configuration. We are 
thus led to consider a wavefront starting at some point 
Po, whose corresponding point po on the cut disk lies 
at a distance tq from the disk center, which maps to 
the cone apex. Experimentally, this scenario can be cre- 
ated deliberately with a colloidal particle that induces 
homeotropic anchoring for the molecules of the liquid 
crystal. We can always cut the cone so that L is exactly 
opposed to Po and then choose coordinates such that L 
lies in the XZ plane. By doing so, we have po = (— ro, 0). 
The geodesic "light rays" may then be parametrized by 
x{X) = — ro + Acoscj, y{X) = Asincj; see Fig.js] Whereas 
computing the geodesies on a generic surface is nontriv- 
ial, for a conical substrate we have a simple analytic map- 
ping of straight lines on the cut disk to geodesies on the 
cone. The associated smectic layers are concentric cir- 
cles centered at po. Note that we could choose the cut L 
along any direction we like as long as we identify the two 
edges. Were we to do so, when a straight line in the Bat- 
tened model hits one of these cuts, we would continue it 
with a straight line emanating from the other cut, mak- 
ing the same angle with the new edge. This ensures that 
the geodesies remain straight and demonstrates that the 
smectic texture is independent of the choice of L. The 
presence of any cusps or grain boundaries in the smec- 
tic does not result from the flattened geometry - all this 
could be computed directly on the cone, for instance. 
Note that we could, alternatively, construct geodesies on 
a (full) two-disk, parameterized by X and Y with induced 
metric: 



ds' = 



cot^ C 



X2 + r2 



dX^ + cot^ C 
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X2 + r2 



dXdY 



1 + cot^ C 



y2 



x2 + r2 



dY' 
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Let X{X),Y{X), Z{X) be the coordinates on the cone 
of the geodesies defined above. A straightforward calcu- 
lation shows that 



where 
k 



r^(A)i 




Y{X) 


= k 







COS ^ X + sin ^ Y 
- sin ^ X + cos ^ F 
cos^Z 



(3) 
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The unit vector field obtained after normalizing this ex- 
pression is given by 



N 



V^l + cot^ C cos2 ^ 



cos(^ - ^) " 
sin(<l> - ^) 
— cos ^ cot ( 



(6) 



Therefore, the projection of N onto the XY plane makes 
an oriented angle <l> — ^ with the X axis (note that ^ 
depends on X and Y through A and uj (Fig. [3|). The 
corresponding vector in the XY plane points along the 
unit direction 



cos(^ - ^) 
sin(<l> - ^) 



(7) 



FIG. 3: Cut-disk view of the cone, showing a parametrization 
of geodesies emanating from a point discUnation. 

The projected layers can also be directly obtained in 
these coordinates by drawing lines which are everywhere 
perpendicular to A/p, with respect to the induced cone 
metric ([2|. Thus, a single prescribed defect, together 
with the constraint of equal layer spacing, uniquely de- 
termines the layer structure everywhere. 

Alternatively, we can compute the layers as the level 
sets of a function D{P) that measures the distance from a 
given point P to the wavefront source Pq. If the cut line L 
is appropriately chosen (so that it contains Pq, for exam- 
ple), the distance between two points P = (X, y, Z) and 
Pq = (Xo, Yq^Zq) on the cone is just the planar distance 
between their counterparts p = {x^y) and po = (xo,^o) 
on the cut disk. In terms of the coordinates ^ (see 
Eq. ([l])), this yields 

D{P) = V(^ + ro)2 + ^2 (8) 
= ^R'^ +rg + 2roi?cos + <l>sinC^. 

The layer structure on the cone, obtained from equally 
spaced level sets of is shown in Fig. |4] The projected 
layers and geodesies seen from above are shown in the left 
panels of Fig [5] Notice that the positive X axis develops 
a grain boundary for any nontrivial deficit angle. 

Indeed, the existence of a cusp is a consequence of the 
Gauss-Bonnet Theorem [8 ; the maximum cusp angle of 
the layers at the grain boundary is equal to n — S and 
occurs at the cone apex. To see this, consider a closed 
path consisting of the geodesic ai uj = ujq from Pq to an 
arbitrary point Pi on the grain boundary, followed by the 
"mirror" geodesic at = —ooq from Pi back to Pq. The 
geodesies form an interior angle ir — ac on the cone at Pi , 
where ac is the cusp angle formed by the layer at Pi . At 
Po, 2ujo is the interior angle formed by the geodesies on 
the cone. Because the geodesic curvature is zero on this 
path, and the integrated Gaussian curvature is simply (5, 
the Gauss-Bonnet Theorem implies 

S = 27T — {tt — ac) — (tt — 2ujo) = o^c + 26^0 (9) 



FIG. 4: Layer structure on the 3D cone for a deficit angle 
S — 57r/4. The image on the right shows the back of the 
image on the left and vice versa. 



(notice that this equation also follows at once from the 
flattened model by elementary geometry). The cusp an- 
gle ac is therefore maximized when we take uq 0+, 
which corresponds to taking Pi arbitrarily close to the 
cone apex. For values of S greater than tt, the cusp angle 
becomes tt at some point to the right of the cone apex, 
meaning that the layers turn back toward the apex. Con- 
sequently, when 5 > 7T the grain boundary is interrupted 
by two new point defects: a +l-index disclination at the 
apex and a — 1-index disclination on the positive X-axis! 
This is shown in the left panels of Fig. |5] and also in 
Fig. [4j The negative-index defect results from the fact 
that the normal of some layers turns through an angle 
greater than or equal to tt, and the outermost such layer 
has a self-intersection on the grain boundary, resulting 
in a locally hyperbolic configuration. By setting F = 0, 
Np = (0, ±1) and solving for X, we discover that the — 1- 
index disclination is located at X = —^^^ro cos ((5/2) 
when this quantity is positive. As (5 ^ 7r+, the —1-index 
disclination coincides with the +1 disclination at the cone 
apex, and for smaller values of 6 the grain boundary is 
free of point disclinations. Note that this disclination 
dipole does not create a dislocation and is an example 
of a pincement lT2\ [T5] that is so "large" as to have gen- 
erated extra internal concentric layers, the dual to large 
Burgers vector dislocations [15^ . 



III. SCHLIEREN TEXTURES 

How would these layer structures appear in an exper- 
iment? In examining nematic and smectic liquid crys- 
talline textures, it is common to view the sample be- 
tween a pair of perpendicularly crossed polarizers. The 
resulting schlieren texture, characterized by dark brushes 
on a bright background, reveals where in the sample the 
molecular orientation aligns on average with the direction 
of either polarizer. For a smectic-A liquid crystal on a 
conical substrate, the molecules are normal to the layers 
aligned along the unit vector field N = {Nx-,Ny-,Nz) 
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FIG. 5: The left panels show the layer structure (black lines) 
superimposed on geodesies (gray lines on the bottom half) 
on the cone, as seen from above (i.e., projected on the XY 
plane), for deficit angles S = 7r/2,7r, and 37r/2, respectively. 
The right panels show the corresponding schlieren textures 
for the same deficit angles. A and Pq label the cone apex and 
the disclination location, respectively. 



in the three-dimensional ambient space. If we were to 
view the sample between a pair of crossed polarizers par- 
allel to the XY plane, we would measure A/^, the nor- 
malized horizontal projection of N. Denoting O as the 
oriented angle between the axis of one of the polarizers 
and the X axis, it follows from Eq. ^ that the intensity 
of the light observed at a point (X, Y) is proportional to 
sin^ [2(^ — ^ — ©)]. Fig. p\ shows the schlieren texture 
along with the layer structure for several choices of the 
deficit angle when = 0. It is interesting to note that, 
besides the defect at Pq (wavefront source), the schlieren 
texture also displays what is usually the signature of de- 
fects, the termination of dark brushes, at the apex and 
at another point farther down the cone at positive X. 
This occurs due to the grain boundary even for (5 < tt. 



FIG. 6: Schlieren textures for deficit angle S = 37r/2, with 
the polarizer direction at angles of G = 0, tt/S, 7r/4 and Stt/S 
with the X-axis, respectively. The analyzer direction rotates 
to remain perpendicular to the polarizer direction. A and Pq 
label the cone apex and the disclination location, respectively. 



when the positive X axis contains no topological defects. 
More surprising deviations from the usual rules governing 
schlieren textures are apparent when J = 37r/2 and we ro- 
tate the polarizers, as shown in Fig.[6j To the right of the 
apex, dark brushes abruptly disappear into the horizon- 
tal axis from below, while other dark brushes spring into 
existence in the upper half-plane, as the polarizers turn 
counterclockwise. In an experiment, such a schlieren tex- 
ture would be the clearest evidence of a grain boundary, 
demonstrating the range of "missing" angles associated 
with a discontinuity in layer normals. Furthermore, the 
number of dark brushes emerging from the point defect 
at the apex is not constant and is odd for certain polar- 
izer angles. In contrast, liquid crystalline textures that 
are continuous except (only) at point disclinations typ- 
ically exhibit a constant, even number of dark brushes 
emanating from each disclination [16 . This strange be- 
havior can be understood by noting that the normalized 
horizontal projection Np of N is not orthogonal to the 
projection of the layers on the XY plane (as opposed to 
the 3D vector N and the layers on the 3D cone which 
are, of course, orthogonal to each other). This can be 
easily seen in the left panels of Fig. [5] and follows from 
the form of ds'^ in ([2). We will come back to this point 
in the next section when we discuss the conical bump. 

A slight generalization of the conical surface above is 
given by a tent, as shown in Fig [7| for S = tt. When 
(5 > TT, a ±l-index disclination pair appears as on the 
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cone, with the +l-index dischnation located at the right 
endpoint of the tent ridge. The layer structure can be 
obtained identically as before by employing the flattened 
model shown in the first panel of Fig ^ 
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FIG. 7: Substrate in the shape of a tent. The top panel shows 
a flattened model from which a tent can be constructed by 
gluing along the lines AB and the lines BC. It also shows a 
point dischnation at Pq and its associated geodesies (dashed 
lines) and layers (solid lines). The second and third panels 
show the corresponding layers in the 3D tent and their 2D 
projection for the case when 6 — tv and AB = APq — 1. A, B, 
and Pq label the cones apices and the prescribed dischnation 
location, respectively. 



IV. SMECTICS AROUND EDGES 

The infinite cone has an isolated singularity at the 
apex. We have also considered the case of a tent, where 
the substrate has an edge, i.e. a line where the surface 
is not smooth. Other examples are shown in the last 
three images of Fig. [l] for which a flattened model is not 
easily obtained because the edge is not straight in the 
xy-coordinate system. The geodesies of such surfaces will 
generally appear kinked at the edge. In fact, an argument 
similar to what is used in geometric optics, in connection 
to Fermat's principle, shows that a geodesic should cross 
an edge following Snell's law. This can be easily seen by 
noting that a geodesic is a curve with constant velocity 
that provides the path of minimal length — and there- 
fore minimal time — between two given points. Note 
that the smectic analog of time is the number of layers 
through which the geodesic passes over a given distance. 
Since the smectic layer spacing is the same on both sides 
of the interface, Snell's law implies that the angle of in- 
cidence equals the angle of "refraction" from the edge, 
where these angles are measured in the tangent planes 
on either side of the edge. The angle of refraction might 
differ from the angle of incidence if an interface separated 
two smectic phases of different layer spacing, as larger 



layer spacing is analogous to smaller index of refraction. 
This could occur in systems of immiscible smectics or at 
first-order transitions between different smectic phases of 
the same material. But we digress. 

Consider the "crater" on the upper right panel of 
Fig. [l] By symmetry, its layer structure (provided some 
boundary condition) can be immediately obtained from 
that on the single cone by reflection across an appropri- 
ate horizontal plane. A more interesting configuration 
is the mountain pass shown in the bottom left panel of 
Fig.[l] Here we can also use symmetry to simplify mat- 
ters. Since the substrate has mirror symmetry across the 
vertical plane that contains the intersection, a geodesic 
that crosses the interface is simply the mirror image of 
a geodesic reflected through the edge. In the flattened 
model of the cone, the intersection will thus appear as 
a boundary F that acts like a mirror, reflecting the "in- 
cident rays" according to the "angle of incidence equals 
angle of reflection" rule on the plane. This is displayed in 
Fig. |8j which also shows the resulting layer structure for 
this case. Notice that a grain boundary is formed at the 
intersection between the surface and the XZ plane for 
all points on the left of the rightmost apex, even between 
the cones. 




FIG. 8: Smectic layers on two intersecting cones with deficit 
angle tt and apices A, B. The layer structure is determined 
by a point dischnation at Pq , located a distance ro from apex 
A as measured on the cone. The horizontal distance between 
the two apices is chosen to be 2ro. The top left panel shows 
geodesies in the flattened out model of the rightmost cone. 
Notice that the cut line L (Fig. [2]) coincides, in this case, with 
the positive and negative y axes so that (0,ro) and (0, — ro) 
represent the same point. The intersection line is represented 
by r. Geodesies on the rightmost cone are represented by 
solid lines while the mirror reflection of those geodesies that 
enter the leftmost cone are depicted by dashed lines. The 
remaining panels show the layer structure on the cones. 

Our discussion so far illustrates the general principle 
that, whenever curvature is present, the constraint of 
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having equally spaced layers leads to singularities in their 
structure, with the appearance of cusps and grain bound- 
aries. In particular, smectics on substrates composed of 
Gaussian bumps have been shown to provide an acces- 
sible system where these ideas take place [8 . We now 
analyze a minimalist and localized version of the Gaus- 
sian bump, the conical bump on the bottom right panel 
of Fig.[ll 
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FIG. 9: Conical bump with deficit angle given by tt. We 
consider as boundary conditions layers parallel to the X axis 
at y ^ — oo. The top panels show the layers on the 3D bump. 
The bottom left panel shows the projected layers (black) and 
projected geodesies (gray) on the substrate; notice that they 
are not orthogonal to each other in the plane metric. The 
bottom right panel shows the projected layers (black) along 
with fictitious layers (thin blue lines) which are constructed 
by demanding orthogonality with the projected geodesies (see 
text). The intersection between the cone and the plane is 
depicted by a dashed circle and A labels the cone apex. 

We consider a situation with the boundary condition 
chosen to be layers parallel to the X axis at Y — oo. 
Before meeting the cone, a geodesic 7 which is normal 
to the layers is a straight line parallel to the Y axis. 
At the interface it deflects according to Snell's law and 
then becomes a geodesic on the cone. There are two 
possibilities at this point, as the geodesies can become 
trapped in the cone or can escape. If 7 enters the cone 
near the Y axis, i.e., with X = Xq close to 0, it will 
reach the X = plane before leaving the cone. As Xq 
grows, 7 will leave the cone and become a straight line 
again before crossing the X = plane. By symmetry, the 
same will happen to the geodesic corresponding to — Xq 
and, as a result, a grain boundary will develop at the 



points of the cone located along the positive Y axis. An 
interesting observation can be made about the escaped 
geodesies. Since a cone is an axisymmetric surface, 7 
must satisfy Clairaut's relation. This means that if p{s) 
is the radial distance of the point 7(5) (in the XY plane) 
from the cone apex and a{s) is the angle that 7^(5) makes 
with a longitude line of the surface, then p{s) sin a{s) 
is constant for each such geodesic. At the boundary of 
the cone all the values of p are the same and a is just 
the angle that 7 makes with the interface. Therefore, 7 
enters and leaves the cone making the same angle with 
the edge. As before, to determine the trajectory of the 
geodesies inside the cone (and their exit point), we may 
either trace straight lines in the flattened model or solve 
the geodesic equations on the cone. This is all illustrated 
in Fig. [9] for a conical bump with deficit angle S = tt. 

We saw in the previous section that the projected lay- 
ers of a non-planar surface are generally not orthogo- 
nal to their projected geodesies. This leads to schlieren 
textures exhibiting an odd behaviour when the sample 
is analyzed with crossed polarizers parallel to the XY 
plane. Now, suppose that we do not initially know that 
the sample is really a curved surface and try to interpret 
it as a planar substrate. Apart from the fact that we 
would be surprised by the unusual pattern of brushes, 
we would also be led to identify a fictitious set of planar 
layers which are everywhere orthogonal to the projected 
geodesies, the latter being inferred from measurements 
under crossed polarizers. Since the projected geodesies 
are not geodesies on the plane by themselves, these ficti- 
tious layers cannot be equally spaced! This is illustrated 
by the bottom panels of Fig.[9| Notice that, as expected, 
both the projected and the fictitious layers agree in the 
planar region but, inside the bump, the latter are highly 
compressed and have the opposite sign of curvature than 
the projected layers. This explains why the schlieren tex- 
tures resulting from a non planar surface look so odd. If 
these fictitious layers were real, they would correspond 
to a high energy configuration, due to compression, and 
therefore would not represent the ground state. The tran- 
sition to an equally spaced structure amounts precisely to 
escaping to the third dimension and assuming the layer 
configuration and shape of the bump. The infinite strain 
in the fictitious layers near the cone apex signals this in- 
compatibility as well. Whether the fictitious layers and 
their geometry can be used as a surrogate to calculate 
the back reaction forces on a deformable surface is an 
open question. 
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